Bull. Inst. Zool., Academia Sinica 27(1): 27-35 (1988) 


MITOCHONDRIAL DNA POLYMORPHISM IN DROSOPHILA 
-~ ALBOMICANS IN TAIWAN! 


HWEI-YU CHANG*?, SU-HUA CHANG**, MEI-LIEN CHU* 
and FEI-JANN LIN? 


* Department of Plant Pathology and Entomology, 
National Taiwan University, Taipei, Taiwan 10765, 
Republic of China 
** Institute of Zoology, Academia Sinica, Taipei, 
Taiwan 11529, Republic of China 


(Accepted November 3, 1987) 
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chondrial DNA polymorphism in Drosophila albomicans in Taiwan. 


Bull. Inst. Zool., 


Academia Sinica 27(1): 27-35. Drosophila alborsicans is one of the most abundant and 
widespread species in Taiwan. The well studied background information makes it 
a good candidate for studying population genetics and evolution. We have used two 
restriction endonucleases, EcoRI and HindIII, to analyze the mitochondrial DNA 
(mtDNA) of 359 isofemale stocks from 29 localities. Those restriction patterns 
reveal that mitochondrial DNA is polymorphic in local populations, but the number 
of major morphs is limited. There are two major EcoRI patterns, and two major 
HindIII patterns. But there are only two major composite EcoRI-HindllI patterns 
instead of four. Different morphs may have different selection values, in other 
words the variation of restriction pattern is probably not neutral. Although the 
mitochondrial DNA is polymorphic within species, the hypothetical primitive 
morph may still be useful information source for reconstruction phylogeny among 


closely related species: 
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Im the past two decades, molecular data 
are frequently used to reveal relationships 
among populations. Compairing to the in- 
formation acquired from proteins, informa- 
tion from DNA is more direct estimates of 
genetic variation. Animal mitochondrial DNA, 
which evolves five to ten times more rapidly 
than single copy nuclear. DNA, in. vertebrates 
(Brown and Wright, 1979; Brown, et al., 
1982; Vawter and Brown, 1986), is a handy 
tool to study relationships among populations 
as well as that among closely related species 
(Avise, et al., 1979a, 1979b, 1986; Ferris, et 
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al. 1981, 1983; Shah and Langley, 1979; 
Yonekawa, et al., 1981; Templeton, 1983). 
Drosophila albomicans is a member of D. 
nasuta subgroup which belongs to D. immigrans 
species group. D. nasuta consists of about 
ten closely related species. Before we can 
use mtDNA data to study the phylogeny of 
those species, we have to know whether there 
is intraspecific variation. The purposes of 
the study are to survey restriction patterns 
of the mitochondrial DNA of the Drosophila 
albomicans in order to reveal the intra- and 
inter-population variation, and to find out 
whether a collection of isofemale stocks 
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maintained in the laboratory can represent 
the natual population. 


MATERIALS AND METHODS 


Flies 


One hundred and fifty three isofemale 
stocks from 29 localities (Fig. 1) in Taiwan 
were maintained in the laboratory for about 
ten years, and two hundred and six isofemale 
stocks from 15 localities were newly estab- 
lished during the year 1985. 


Extraction of mitochondrial DNA 


Mitochondrial DNA was extracted from 
larvae according to the protocol modified ` 
from that reported by Coen, ef al. (1982). 
Briefly, 15 larvae were gently homogenated 
in 320 ul of solution I (10 mM Tris, pH 7.8, 
60 mM NaCl, 5% sucrose, and 10mM EDTA), 
and lysed by the addition of 400,41 of 
freshly mixed solution II (300 mM Tris, pH 
9.0, 1.25% SDS, 5% sucrose, 100 mM EDTA, 
and 8 ul/ml diethylpyrocarbonate) After the 
reaction mixtues were incubated at 65°C for 


Fig. 1. Twenty nine collecting sites of Drosophila albomicans are indicated on the 


map of Taiwan. 


See Table 3 for the names of these localities. 
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30 min, 120 wl of 3 M sodium acetate was added. 
Stand on ice for 45 min, and then centrifuge 
on Eppendorf for 10 min at 4°C. DNA was 
precipitated by isopropanol at room tempera- 
ture, and resuspended in 250 wl double dis- 
tilled water containing 0.6 ul diethylpyrocar- 
bonate (freshly mixed). Stand at room tem- 
perature for 20min, and then add another 
250 u] double distilled water. Mitochondrial 
DNA. was precipitated by ethanol and dried 
under vaccum. 


Digestion of restriction endonucleases and 
agarose gel electrophoresis 


Mitochondrial DNA, extracted by the 
method described above, was dissolved in 20 
ul double distilled water. The digestion of 
restriction endonucleases was performed ac- 
cording to the procedure described by Maniatis, 
et al. (1982). The reaction was stopped by 
the addition of 1/10 volumn STOP solution 
(0.05% bromophenol blue, 50mM EDTA, 
0.05% SDS, 50% glycerol, and 2 mg/ml boiled 
RNase). The restriction pattern was analyzed 
on a 0.8% agarose gel in Tris-acetate buffer 
(4 mM Tris, pH 8.1, 2mM acetic acid, and 


Hind IIl 


0.2mM EDTA). After staining with ethi- 
dium bromide, DNA bands were visualized 
under UV light. Size of fragments were 
estimated by a graphical method using HindIII 
digested lambda phage DNA as size standard. 


RESULTS 


MtDNA from 359 isofemale stocks was 
subjected for HindIII and EcoRI digestion. 
Since the enzyme digestion was done by 
using crude extract of mtDNA, individual 
results on gel are not good representatives 
of restriction patterns. Besides, the DNA 
fragments ranging from 11.8 kb to 1.1 kb do 
not appear clearly on a same gel (Fig. 2). 
We use the average of about 30 measure- 
ment data for each restriction pattern to 
produce the graph of Fig. 3. Size determina- 
tion was performed by running mtDNA 
restriction fragments, mtDNA restriction 
fragments+DNA size standard, and DNA size 
standard side-by-side on a gel. Size of large 
fragments was determined on low concentra- 
tion agarose gel, and that of small fragments on 
high concentration gel. There are 7 restrition 
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Fic. z. Fhoto of seven HindIII restriction patens and 4 EcoRI restriction 
patterns of mitochondrial DNA. 
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Fig. 3. Restriction patterns are showed 


indicated beside. 


patterns (A, 61.0%; B, 35.3%; C, 1.1%; D, 
0.6%; E, 03%; F, 0.3%; and G, 1.4%) 
resulted from the digestion with HindIII, 
and 4 (A, 30.6%; B, 66.1%; C, 2.5%; and D, 
0.8%) with EcoRI (Fig. 3). If the mtDNA 
in D. albomicans evolves from the same 
ancestor, we may construct a network to 
connect these restriction patterns. These res- 
triction patterns have been arranged in Fig. 
4 so as to minimize the number of restric- 
tion sites that need to be gained or lost in 
network. In all cases the transition from 
one pattern to another requires just a single- 
site change (gain or loss). As shown in 
Fig. 4, these two restriction endonucleases, 
HindIII and EcoRI, yield a total of 13 
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graphically with sizes of fragments 


mtDNA restriction sites in the 359 isofemale 
stocks; 4 sites (31%) are shared by all 
stocks, the other 9 (6927) are polymorphic. 
Pattern A of HindIII digestion (219 stocks) 
and pattern B of EcoRI digestion (237 stocks) 
are the main patterns found in D. albomicans 
on this island. 

All of the composite patterns (morphs) 
found are showed in Table 1. Only 14 out 
of 28 hypothetical morphs are observed in 
this study. The frequency distribution of 
those morphs does not fit the expected fre- 
quency calculated from single digestion data 
(Table 1. Two mtDNA morphs, 2 and 4, are 
common. Morph 2 is presented in 205 (57.127) 
and morph 4 in 92 (526%) of the 359 D, 


MITOCHONDRIAL DNA 31 


Hind Ill 


G 


$-G-6-09-0 


EcoRI 


© 


Fig. 4. Diagram of the restriction patterns obtained with endonucleases 
HindIII and EcoRI in D. albomicans mtDNA. The most common 


restriction patterns are HindIII A and EcoRI B. 


Locations 


shown for the restriction sites are arbitrary. 


albomicans isofemale stocks. The remaining 


twelve morphs consist of just 17.3% of the 


population. Six of them are found in single 
stocks, and the other six are scatteredly dis- 
tributed. 

Table 2 shows the distribution of 
mtDNA morphs in old laboratory maintained 


and newly established isofemale stocks. Ac- 
cording to x? test, there is no significant 
difference between these two groups. So the 
data of old stocks and that of new stocks are 
combined and the result of this survey is 
summarized in Table 3. MtDNA seems to 
be polymorphic in local populations, and to 


TABLE 1 
Restriction pattern for the 14 mtDNA morphs found in D. albomicans 
IE 


Restriction battern Stocks Expected 
Morph HindIII EcoRI Number Percent percent 
1 A A 13 3.6 18.67 
2 A B 205 57.1 40.32 
3 A C 1 0.3 1.52 
4 B A 92 25.6 10.80 
5 B B 24 6.7 23.33 
6 B C 8 2.2 0.89 
7 B D 3 0.8 0.28 
8 C A 3 0.8 0.34 
9 C B l 0.3 0.73 
10 D A 1 0.3 0.18 
11 D B EI 0.3 0.40 
12 E A 1 0.3 0.09 
13 F B 1 0.3 0.20 
14 G B 5 1.4 0.92 


RE 和 ES 
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; TABLE 2 | 
The distribution of mtDNA morphs in old laboratory maintained and 
newly established isofemale stocks of D. albomicans 


Morphs 
Stocks 1 2 3 4 5 6 7 8 9 10 1] 12 13 14 
Old 5 90 1 4] 8 1 2 0 0 1 1 i 0 2 
New 8 115 0 51 16 7 1 3 1 0 0 0 i 3 
Total 13 205 1 92 24 8 3 3 1 1 1 l 5 
TABLE 3 


The distribution of mtDNA morphs in D. albomicans 


No. of Morphs 
Locality stock 1 2 3 4 5 6 7 8 9 10 IH 12 13 14 
1. Baileng (G4) 16 11 3 1 1 
2. Chia-yi ( 嘉 闵 ) 27 1. .11 8 3 3 1 
3. Dawu CX) 4 3 1 
4. Dungpu (Xii) 7 5 5 2 
5. Funglin (JEU) 17 10 ' 4 2 1 
6. Guguan (BB) 5 E 4 1 
7. Jrben ( 知 本 ) 3 . 2 1 
8. Juolan (Wf) 17 1 9 4 3 
9. Kending (ÆT) | 42 1 22 12 4 1 2 
10. Lanyu (Hti) E 1 1 l 
1. Majia GEX) "ES 1 | 
12. Meishan (梅山 ) 12 1 6 4 ] 
13. Muguashan (AJM) 6 1 5 
14. Nan-gang (南港 ) 5 3 2 
15. Puli ($E) 17 7 6 3 1 
16. Rueifang (35525) 4 3 1 
17. Rueisuei (瑞穗 ) 18 10 6 1 1 
18. Sandimen (hf) 4 1 2 1 
19. Sanshia (=k) 4 i 3 
20. Shanhutan (FWA) 7 6 Í 
21. Shitou (7) 9 9 
22. Shrtoushan (38) 10 . 1 4 3 2 
23, Taibei (XL) 1 1 
24. Tailuge (KSI) 24 12 9 2 1 
25. Taiping (KA) 1 1 
26. Wufeng (XX) 1 fe’ 
27. Wulai (Esk) 48 27 ] 17 1 1 1 
28. Yilan CES) I 1 


29. Yuenshuel (Jzjk) 46 5 28 8 2 1 1 1 
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have a couple of common major morphs in 
all populations. 


DISCUSSION 


The size of the D. albomicans mtDNA is 
estimated to be 16.5 kilobases (kb), which is 
about average for animals. The range for 
most animals studied is 15.7 to 19.5 kb 
(Brown, 1983). No length polymorphism was 
Observed among the 359 isofemale stocks. In 
the strict sense, no size difference which is 
larger than 0.5 kb has been found in this 
study, but the restriction digestion shows 
polymorphism. The laboratory maintained 
isofemale stocks are still good representa- 
tives to natural population on the mtDNA 
variation. 

Five localities, each represented by only 
one stock, do not provide useful data for 
analysis of local population variation. All 
the remaining populations, execpt the one at 
= Shitou including 4 new stocks and 5 old stocks, 
are shown to be polymorphic. According to 
our observation, even for the flies collected 
from a single bait there usually exists poly- 
morphism. One possible explanation for the 
monomorphism in Shitou is that this local 
population has been maintained in small 
size and is insulated, since the possibility of 
having the same morph in 9 collections 
from 1971 to 1985 is low. Genetic dirft plus 
the inhibition of migration may be the cause 
of monomorphism. 

The fact of the high polymorphism in 
mtDNA as well as the scattered distribution 
of minor morphs in populations is conflicting 
to the common thought that mtDNA is 
homogenous within a small population be- 
cause of its clonal inheritance (Brown, 1983). 
This result may be partly explained by a 
high mutation rate of mtDNA. Therefore, 
the presence of a certain restriction site may 
be the result of convergent or parallel evolu- 
tion instead of common ancestry of popula- 
This phenomenon is quite different 
Because 


tions. 
from chromosomal inversion data. 


of the extremely low mutation rate of varia- — 


ble inversions, a population loses or gains a 


chromosomal inversion through genetic drift 
and migration, respectively. Thus, the presence 
of a particular chromosomal inversion type 
instead of a particular mtDNA restriction 
site may reveal geneology of populations. 

The analysis of nucleotide substitutions 
between pairs of DNA Sequences can not 
directly reveal the phylogenetic relationship of 
species simply because of polymorphism. This 
problem may be solved by increasing sample 
size of each population and by using a large 
number of restriction enzymes as recom- 
mended by Nei and Taijma (1983). However, 
the phylogenetic tree of populations cons- 
tructed from the frequency distribution data 
of restriction patterns may still be screwed 
up by selection and nonrandom mutation, 
even if the sample size is large enough. By 
using a large number of restriction enzymes 
may help us minimize the bias resulted from 
selection and nonrandom mutation. 

The other alternative is that the phylo- 
genetic relationship of species may be recons- 
tructed by the primitive morphs of mtDNA. 
How to recognize the primitive morph will 
be a difficulty in this aspect. In our case, 
we consider HindIII pattern B and EcoRI 
pattern B to be the primitive forms (Fig. 4), 
then morph 5 (Table 1) should be the pri- 
mitive composite form. If this is true, the 
most abundant form is not the most primi- 
tive form. The two major morphs in the 
present populations are both derived from 
this hypothetical primitive form. This result 
may also suggest that the hypothetically pri- 
mitive form is not the fittest on this island. 
The fact that all of those large samples have 
the same major morph in spite of the 
population differentiation showed by our 
previous study (Lin and Chang, 1986) in- 
dicates this major morph really have advan- 
tage on fitness in Taiwan. 8 

In this study, there are six morphs 
which appear in single stocks. The other six 


minor morphs ranging from three to 24 
stocks show scattered distribution. In D. sub- 


obscura, however, six minor morphs appear 
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to be narrow endemic, each found in only 
one locality (Latorre, et aL, 1986). This 
inconsistancy may be due to sample size. 
Five out of the six minor D. subobscura 
morphs were found in single strains. And 
the only informative morph appears in just 
two strains. 

In conclusion: (1) mtDNA is polymor- 
phic within local populations; (2) mtDNA 
does not show strong differentiation among 
populations; (3) there are two major morphs 
presented in D. albomicans on this island; 
and (4) neither of the two major morphs is 
the hypothetical primitive morph. 
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Xr Æ E (Drosophila albomicans) WB 18 Jk 2 
< Bit WR BS DNA Z BBE 


RHE (Drosophila albomicans) 在 训 湾 数量 多 且 分 侯 广 ， 而 且 已 有 许多 研究 资料 ， 因 此 很 通 宜 用 
来 研究 族 硬 遗传 与 进化 。 我 们 探 用 两 种 限制 酶 EcoRI 和 HindIII ROHR BS 29 个 地 区 的 359 个 
ESHER du RES DNA o f REURBI EEA» ARES DNA 仍然 呈 多 驴 性 。 EcoRI 限制 
酶 切 出 两 种 主要 的 型 式 ，HindIJII gh > m EcoRI-HindIII 的 组 合 限制 酶 切 型 却 只 有 两 种 主要 的 型 
A> HMDA SERRA LWER > Bez? RARE aH EP ER BER 
种 内 粒 腺 体 DNA 具 多 能 性 ， 但 假设 的 原始 限制 酶 切 型 对 於 近 各 种 的 种 问 关 保 之 建立 仍 可 提供 有 用 的 资 
to 


